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SUMMARY 


A reevaluation of existing flight data obtained for a supersonic 
canard missile configuration of the boost-glide type has been made to 
determine derivatives not previously evaluated. These derivatives, 
together with those previously evaluated and published, have been 
utilized to determine some typical airframe frequency responses based 
on the three -degree -of -freedom longitudinal equations of motion. 

For constant flight conditions, it is shown that besides the usual 
drag-coefficient variation, the velocity derivatives also vary with 
angle of attack or Mach niomber. The effect of angle-of -attack variation 
on the missile frequency responses which include the velocity derivatives 
in their solution is appreciable only in the low-frequency region of 
operation (below 10 radians/sec ) . V/hen the velocity derivatives are 
neglected, this low-frequency variation with angle of attack is different. 
Some of these differences are pointed out in the results to indicate the 
significance of including the velocity derivatives in the solution for 
the transfer functions. 


INTRODUCTION 


Increased interest in the phugoid motion has been expressed from 
various sources, especially in connection with automatic airspeed or 
altitude control systems. The present study of the three-degree-of- 
freedom motion actually arose out of an investigation of an automatic 
altitude control system wherein the altitude- sensing instruments operate 
at a frequency considerably below the high-frequency airframe mode. In 
this situation, the influence of the low-frequency airframe dynamics on 
the overall stability of the system may prove to be important. 
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The main purpose of the present paper is to evaluate the longitudi- 
nal transfer functions of a supersonic canard missile by means of a 
reevaluation of the data obtained from previous flight tests of the mis- 
sile in order to include the variation of drag coefficient with angle of 
attack and the variation of force and moment coefficients with velocity 
in the solution of the three-degree-of -freedom equations of motion. The 
variation of these additional derivatives with Mach number is presented 
herein for several values of angle of attack. The longitudinal transfer 
functions in the form of frequency-response plots for several flight 
conditions are also given. 

In reference 1, the effect of a fourth degree of freedom (namely, 
variation of height) on the longitudinal motion of full-scale aircraft 
is discussed. Although this additional degree of freedom may influence 
the longitudinal motion of missiles to some extent, the study of this 
effect was not included in the present investigation. The present inves- 
tigation is mainly concerned with the additional effects of variation of 
force and moment coefficients with velocity, as compared with previous 
three -degree -of -freedom studies (ref. 2, for example) wherein these 
effects were not considered. 


SYMBOLS 


9 pitch-attitude angle, radians 

a angle of attack, radians unless otherwise indicated 

P angle of sideslip, radians 

7 flight -path angle, radians 

^ critical damping ratio 

5 canard elevator deflection, radians 

m mass, 5*05 slugs 

W weight, 165 pounds 

S wing area, 4-.1 sq ft 

Y stability axis which passes through center of gravity and 

is perpendicular to vertical plane of symmetry 

ly moment of inertia about Y-axis, 51-5 slug-ft^ 
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c 

^sm 


M 

V 


u 


V 


mean aerodynamic chord, 1.776 ft 
static margin, ft 

Mach number 
velocity, ft /sec 


q 

g 

t 

D 

KG(D) 

AR 

db 

PA 

CD 

Cl 

Cd 


dynamic pressure^ lb /sq ft 
acceleration due to gravity^ ft /sec^ 
time^ sec 

differential operator^ d /dt 

transfer function expressed as a linear fiinction of D 
amplitude ratio of KG(D) 
decibels, 20 log^QAR 

phase angle of KG(D), deg 
angular frequency, radians /sec 
damped natural frequency, radians/sec 
lift coefficient, Lift /qS 
drag coefficient , Drag/qS 

pitching-moment coefficient , , moment 

qSc 


CL5 = 

Ci^ = 6CL/aa 
Cm5 = ^Cjjj|d6 
CiiIq^ = 



CONFIDENTIAL 



4 


COKFIDENTIAL 


NACA RM L54C02 


Cm = SCjh 

^ 2V 

CDq^ = hC-Q^hcL 

CDu = 

^Lu = 

Cmu = 

Subscripts : 
p phugoid 

n natural frequency 

o equilibrixun value 

Dpt over a symbol denotes derivative with respect to time. 

METHODS AND APPARATUS 
Equations of Motion 

The three-degree-of -freedom longitudinal equations of motion, 
assuming small disturbances from steady symmetric horizontal flight, are: 

(is ’’o)® * ^ ^ + 2CDju . 0 

^ =in rc)e - ^ - (=iu + = ° 
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Various derivations and forms of these equations can be found elsewhere 
(see, for example, refs. 2 to 7). In these references various notations 
are used; in some more terms are included and in others terms are neg- 
lected which are included here. The equations presented here, however, 
satisfactorily define the motion of the missile under consideration with 
the assumption of small disturbances. The notation used is based on the 
stability axis system and parallels the more widely used two-degree -of - 
freedom equations (see ref. 8). 


Flight -Data Presentation 

Airframe .- The airframe on which this analysis is based is an all- 
metal research model of the canard missile type. A photograph and a 
sketch of the actual model used for previous flight testing is shown in 
figure 1. References 9 and 10 give more complete descriptions of this 
model and present the flight data in the form of aerodynamic stability 
derivatives used in previous two-degree-of -freedom analyses. In refer- 
ence 8, values of these derivatives are tabulated for typical flight 
conditions. The additional derivatives necessary for the three-degree- 
of -freedom study presented herein were obtained from a reevaluation of 
the flight data presented in references 9 and 10 and from the results of 
other flight testing not previously published. • 

CDa*- This derivative. was evaluated from the results of flight tests 

not previously reported, wherein the model was roll stabilized and pulsed 
in pitch and yaw, and the values of Cjj based on total drag along the 

velocity vector were plotted against total angle of attack \/a2 + p2. 
Sufficient curves were obtained to cover a Mach number range of approxi- 
mately 1 to 1.7* The slopes of tangents to curves obtained in the fore- 
going manner were then measured to determine the Cp^ data presented 

herein. Since a symmetrical cruciform missile is being dealt with, a 
value of Cp based on total drag and total angle of attack is equiva- 

lent to one based on measurements in the pitch plane alone and is valid 
for use in the longitudinal equations of motion. 

^Du^ ^Lu^ this case the variations of the coeffi- 

cients Cp, Cp, eind Cm with velocity were first obtained. Sufficient 
data were available to cover a velocity range of approximately 1,000 
to 2,000 feet per second, assuming standard sea-level conditions. The 
angle-of -attack range covered was as high as 10° or 12°. The slopes of 
tangents to these curves were then measured to obtain the nondimensional 
derivatives such as Cp^ = 8Cp/8u which are presented herein as functions 

of Mach number for constant values of angle of attack. 
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Frequency Responses 

A solution of the equations of motion to obtain the transfer func- 
tions can only be obtained by assuming an equilibrium Mach number and 
then selecting values of and 0^^ for the trim angle- 

of -attack value under consideration. Solutions based on the foregoing 
assumptions are presented in the results as frequency-response plots of 
the transfer functions. A particular transfer function, for instance 
e/6, is obtained by expansion of its determinant to yield the function 


I = KG(D) 


The frequency response of this function is then plotted through the use 
of templates as described in chapter 8 of reference 11. 


RESULTS AND DISCUSSION 
Derivatives 


Cp .- Plots of Cp„ against Mach number for constant values of a 

are shown in figure 2. For a constant Mach number, this derivative exhib- 
its a decided increase with increasing a and for the lower values of a, 
Cpa decreases to a steady value with increasing M. The decrease with 

increasing M is more apparent at the higher values of a for the Mach 
number range shown. 

CPu, Cpu, and Cnl^.- The velocity derivatives Cp^, CLu, and Cm^ 

are plotted against M ^for constant values of a in figure 5- Tbe 
trends exhibited by Cp^ and Cl^ (figs. 5(a) and (b)) are somewhat 

similar in that at subsonic Mach numbers the curves show positive values 
which decrease to zero at approximately M = 1.05 and then increase 
negatively to maximum values between M = 1.2 and M = 1.55* From this 
point, a more gradual rise back toward zero is noted for Cp^ than for 



Curves of against M for constant values of a are shown in 

figure 5(c). From these curves it is seen that C^^ is negative at sub- 
sonic Mach numbers, and is zero for all values of a at approximately 
M = 1.05. The values of Cm^_^ then increase to maximum positive values 
at M = 1.55 approximately, after which a decrease with increasing M 
is shown. The uniform spread in the curves of this figure indicates that 
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the variation of this derivative with a for a constant Mach number is 
fairly linear. 

Cj) and Cl-- The values of Cj) and Cl used in confuting the 

transfer functions presented herein are shown in table I. Except for 
the zero-lift drag coefficient^ the Cj^ values have not been previously 

published. The Cl values are presented for completeness although they 
are indirectly available in published form elsewhere. 


Frequency Responses 

Typical missile frequency responses obtained from a solution of the 
three-degree-of-freedom equatiohs of motion are shown in figures it- and 5. 
These responses are based on standard sea-level conditions at M = 1.6 
with a missile static margin of 0.29kc. The responses 0/6, a/b, 
and u/6 were obtained directly from a solution of the equations of 
motion as given herein previously, while the substitution 7 = 0 - a 
was employed to obtain the 7/6 responses (figs. i|-(d) and 5(d)). For 
each of the variables, the amplitude and phase response obtained for 
trim angle-of -attack values of 0°, 4°, and 10° are shown. This implies 
that the derivatives and Cm^j^ and the values of Cp 

and Cp correspond to the trim angle- of -attack value under consideration 
for each solution, despite the variations with cx shown in the previous 
figures and table I. 

Effect of including velocity derivatives .- The frequency responses 
obtained at. M = 1#6 including tbe etifeetp of* velocity derivatives are 
shown in figure 4. Phase shifts as the frequency approaches zero and 
changes in static sensitivity are shown for the trim angle-of -attack 
values investigated. In figures 4(a) and 4(d), for example, as the fre- 
quency decreases the 0/6 and 7/6 phase c\irves are shown to tend 
towards -90°^ 0 , and I80 as the trim angle of attack is changed from 
0° to 4° to 10°, respectively. At zero trim a the solution is similar 
to that obtained for two degrees of freedom with an integration effec- 
tively present due to the zero lift force. At a = 4° the character- 
istic equation (denominator of transfer function) is of the biquadratic 
form 


(d2 - 4 - 2^pO)^D + -f 2^cDnD + ay/j 

with a positive static sensitivity yielding a zero initial phase relation. 
For a = 10°, however, the numerators of the 0/6 and 7/6 transfer 
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functions each contain a negative linear term yielding a negative static 
sensitivity and causing the 180° initial phase shift. 

For angular frequencies above 10 radians/sec, with the exception of 
u/5 (fig. ^(c)), changing the trim angle of attack has only a slight 
effect on the frequency response obtained, since for this particular 
missile the static stability derivatives have previously been shown to 
be fairly linear. The effect of changing the trim angle of attack would 
therefore be considered negligible in control- system applications where 
the controlling elements are designed to operate at a frequency in excess 
of the high-frequency airframe mode. The reasoning in this case is that 
the overall transient behavior of a control system in combination with 
the airframe woiuLd be essentially the same regardless of the trim a, 
since the contribution of the airframe would not change appreciably with 
angle of attack in the region of the control- system operating frequency. 

The curves presented in figure 4 also indicate that the phugoid mode 
is more heavily damped than the high-frequency mode. The values of 
phugoid critical damping ratio given in the first part of table II are 
0.84 at a = 4° and I.96 at a = 10°, compared with values of ^ O.I5 
for the high-frequency mode. The large value of at a = 10° indi- 

cates that the phugoid mode in this case is nonoscillatory, actually con- 
sisting of two aperiodic modes. This large value of can be attrib- 

uted to the large drag obtained for supersonic flight at this angle of 
attack. The values of damped natural frequency given in table II where 
the velocity derivatives are included are approximately 0.15 radian/sec 
for the low-frequency mode and 26 radians/sec for the high-frequency 
mode. The values of the static sensitivity (actual static ratios) of 
the various transfer fxmctions are also given in table II. With the 
exception of a/5, the static sensitivities are infinite at a = 0° and 
finite when a value is assumed for a. In any case the absolute values 
of the 0/6 and 7/5 static sensitivities are fairly large (for 
instance, the static value of 0/5 is 29-l at a = 4° and -l4.8 at 
a = 10° and the infinite value at a = 0° is, of course, due to the 
previously mentioned integration.) These large or negative values of 
static sensitivity, however, may not be too significant since, with the 
assumption of small perturbations, large steady-state values would o nl y 
be applicable for extremely small values of 5. 

Effect of neglecting velocity derivatives .- The frequency responses 
shown in figure 5 are comparable to those shown in figure 4 except that 
in figure 5 the velocity derivatives are neglected (that is, 

= C];^ = Cjji^ = 0^. The characteristic-equation and static-sensitivity 

values for this case are also summarized in table II. Except for varia- 
tions in static sensitivity, table II, however, does not reveal anything 
very significant concerning the effect of the velocity derivatives. 
Referring again to figures 4 and 5> the variations between the responses 
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shown in figure 5 and the corresponding responses shown in figure 4 are 
mainly in the low-frequency region. A close comparison of the 
0/5 responses shown in figures 4(a) and 5(a) reveals some differences 
in the shapes of the amplitude and phase responses in the low-frequency 
region. The nature of the responses, however, makes the effect of 
including or neglecting the velocity derivatives more readily apparent 
in a comparison of the a/5 responses (figs. 4(b) and 5 (d))* In fig- 
ure 5(D) it is shown that changing the angle of attack has very little 
effect on the response. This result arises from the mathematics of the 
transfer function, since the numerator terms practically cancel the 
phugoid quadratic regardless of the angle of attack. However, when the 
velocity derivatives are included, as in figure 4(b), this is not the 
case. Here the amplitude-ratio curves exhibit noticeable change in 
static sensitivity when a value of a is assumed and the phase curves 
show leading characteristics in the low-frequency region. The leading 
characteristics arise from the fact that the numerator terms break at 
lower frequencies than the phugoid mode of motion. 

Effect of Mach mamber change .- The data presented in table III are 
comparable in every respect to those presented in table II except that 
the values presented are based on a Mach number of 1.2. Variations 
between the corresponding values of the two tables can be noted. Perhaps 
the most significant of these is the value of the phugoid critical damping 
ratio at a = 4^ for the case where the velocity derivatives are 
included. In this case = 0.574 and the low value can be attributed 

to a high value of Cd^. Reference to figure 5(a) shows that for the 
a = 40 

curve, Cp)^ -0.07 or very nearly its negative maximum for this 
angle of attack. The absolute value of Cp)^ for this case is actually 
greater than Cj^ (see table l) . Although the curves are not shown, the 
frequency responses for M = 1.2 and a = 4^ would exhibit more peaking 
of the amplitude ratio and more rapid phase changing in the region of the 
phugoid damped natural frequency (0.II5 radian/sec) than any of those 
shown for M = 1.6 because of the lower critical damping ratio. 


Three-Degree-of-Freedom Approach to Future 
Automatic-Control Studies 

Since the frequency response obtained with constant flight condi- 
tions is shown to vary with the trim angle of attack, particularly in the 
low-frequency region, the analysis of automatic -control systems which 
operate at low frequencies — for example, altitude or airspeed controls — 
should define the airframe on the basis of three degrees of freedom. The 
transfer-function approach may be essential in preliminary analysis 
employing linear servomechanism theory. Therefore, in analyzing the 
altitude control, for example, depending on the angle -of -attack range 
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anticipated, three or more airframe transfer functions should be com- 
puted for each case with flight conditions and all other system compo- 
nents held constant. Then, depending on the con 5 )onent gain adjustments, 
it is probable that the stability of the system can be adequately pre- 
dicted if no instabilities arise, regardless of the trim angle of attack 
on which the airframe transfer function happens to be based. It is 
obvious that more than one transient solution will be obtained for ein 
equilibrium Mach number under these conditions. Although it can o nl y be 
surmised at this time, it is anticipated that the actual transient 
behavior of the entire control system will be within the range of the 
solutions obtained for a number of different trim values of angle of 
attack. In order to obtain the unique solution for a constant set of 
flight conditions, the airframe can be represented by equations of motion, 
wherein some of the coefficients are represented as functions of the vari- 
able a. The airframe defined in this manner is nonlinear but can be com- 
bined with the rest of the automatic-control-system components through 
the use of an analog computer (REAC, for example) to obtain the unique 
solution. This, however, is a lengthy process; nevertheless, it may 
eventually be necessary to determine the validity of the previous 
surmise. 


CONCLUSIONS 


The additional derivatives necessary for solution of the three- 
degree-of -freedom longitudinal equations of motion, namely, the velocity 
derivatives and the slope of the curve of drag coefficient against angle 
of attack, are shown to vary with angle of attack. This variation affects 
the missile frequency responses mainly in the low-frequency region. 
Therefore, in the analysis of a control system which operates at a high 
frequency, the contribution of an airframe with linear static character- 
istics would be essentially the same for all angle-of -attack values. 

A comparison of the frequency responses obtained when the velocity 
derivatives are included with those obtained when the velocity deriva- 
tives are neglected reveals variations in the low-frequency region. The 
angle -of -attack frequency response, for example, exhibits noticeable 
changes in static sensitivity and more pronounced leading phase charac- 
teristics when the velocity derivatives are included than when the 
velocity derivatives are neglected. 

The critical dan 5 )ing ratio was found to be higher for the phugoid 
mode of motion than for the high-frequency airframe mode. For low super- 
sonic Mach numbers, however, where the derivative of the drag coefficient 
with respect to velocity is near its negative peak, the phugoid critical 
damping ratio may assume a fairly low value. 
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For the supersonic canard missile used in this investigation, the 
velocity derivatives reverse sign in the transonic region and have their 
maximum absolute values for supersonic flight in the range of Mach number 
from 1.2 to 1.55* Their absolute values then decrease with further 
increase in Mach number. 

In applying three -degree-of -freedom transfer functions to automatic- 
control-system analysis, more than one airframe transfer function should 
be investigated for an equilibrium Mach number and constant automatic- 
control component gains. Then it is probable that the stability of the 
overall system can be adequately predicted if no instabilities arise, 
regardless of the trim angle of attack on which the airframe transfer 
function happens to be based. It is further surmised that the unique 
transient solution for an equilibrium Mach number will be within the 
range of the solutions based on airframe transfer functions covering a 
rajige of assimied trim angle-of -attack values. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February l6, 195^ • 
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TABLE I 


Cd and Cl values 


a, deg 

Cl for — 

Cl for - 

M = 1.2 

M = 1.6 

M = 1.2 

M = 1.6 

0 

0.05 

O.OA 

0 

0 

A 

.065 

.051 

.211 

.182 

1 

10 

.152 

.128 

.527 

.456 
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TABLE II 


SUMMARY OF TRANSFER-FUNCTION CHARACTERISTIC 


EQUATIONS AND STATIC SENSITIVITIES 


[m 


1.6; sea level; 



Characteristic Equation 


Including velocity derivatives 

a, deg 

radians/sec 


cDjj, radians/sec 


0 


„ 

26.25 

0.152 

k 

0.132 

0.84 

26.2 

.152 

10 

L 

.207 

1.96 

26.1 

.145 

Neglecting velocity derivatives 

0 


_ „ _ „ 

26.25 

.152 

4 

.105 

.83 

26.2 

.152 

10 

.166 

1.31 

26.2 

.152 


Static Sensitivities 


Including velocity derivatives 

a, deg 

e /6 

a/6 

u/6, per radian 

7/6 

0 

00 

0.9 

00 

00 

4 

29.1 

.5 

- 4.2 

28.6 

10 

- 14.8 

.5 

-1.7 

-15.5 

Neglecting velocity derivatives 

0 

00 

.9 

00 

00 

4 

^ 41.8 

• 9 

-6.7 

40.9 

10 

-28.9 

.9 

- 2.6 

-29.8 
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TABLE III 


SUMMARY OF TRANSFER FUNCTION CHARACTERISTIC 
EQUATION AND STATIC SENSITIVITIES 

[m = 1.2; sea level; = 0.339c] 


Characteristic Equation 


Including velocity derivatives 

a, deg 

radians/sec 


0^, radians/sec 


0 



22.8 

0.146 

4 

0.115 

0.37^ 

22.8 

.146 

10 

.182 

1.23 

22.7 

.142 

Neglecting velocity derivatives 

0 



22.8 

.14 

4 

.114 

.711 

22.8 

.146 

10 

.179 

1.07 

22.7 

.146 


Static Sensitivities 


Including velocity derivatives 

a, deg 

0/5 

a/8 

u/5, per radian 

7/8 

0 

00 

.8 

00 

00 

4 

8 

.6 

-5.6 

7.4 

10 

-11 

.5 

-2.2 

-11.5 

Neglecting velocity derivatives 

0 

00 

.8 

00 

CO 

4 

26.6 

.8 

-5.7 

25.8 

10 

-16 

.8 

-2.3 

-16.8 
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Figure 1.- Supersonic research model. 
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M 

Figure 2.- Variation of Cd with Mach ntmiber for constant values of 

OG 

ajigle of attack. 
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(t>) 

Figure 3»- Continued. 
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(a) e/6. 

Figure k.- Three-degree-of-freedom longitudinal frequency responses 

Including effect of velocity derivatives for three values of trim a. 
These responses are based on standard sea-level conditions at M = 1.6 
with Xgjjj = 0 . 294 ^. 
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(b) a/6. 

Figure 4.- Continued. 
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cj , racha/)s/sec 


5 /O 


50 too 


(c) -u/6. 

Figure k.- Continued. 
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(d) 7 / 6 . 

Figure 4.- Concluded. 
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cj ^ radians/ sec 


5 !0 


50 too 


(a) 0/6. 

Figure 5»- Three-degree-of -freedom longitudinal frequency responses 

neglecting effect of velocity derivatives for three values of trim a. 
These responses are based on standard sea-level conditions at M = 1.6 
with Xsm = 0.294c. 
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(b) a/6. 

Figure 5*- Continued. 
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(c) -u/6. 

Figure 5*~ Continued. 
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(d) 7/6 . 

Figiire 5 .- Concluded. 
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